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Abstract 

In seeking a model solving the coincidence problem, the effec- 
tive Yang-Mills condensate (YMC) is an alternative candidate 
for dark energy. A study is made for the model up to the 2- 
loop order of quantum corrections. It is found that, like in the 
1-loop model, for generic initial conditions during the radiation 
era, there is always a desired tracking solution, yielding the cur- 
rent status Qa — 0.73 and Q m ~ 0.27. As the time t — ► oo the 
dynamics is a stable attractor. Thus the model naturally solves 
the coincidence problem of dark energy. Moreover, if YMC de- 
cays into matter, its equation of state (EoS) crosses -1 and takes 
w ~ —1.1, as indicated by the recent observations. 
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1. Introduction 



The energy content of the present Universe is such that ~ 0.73 and Q m ~ 0.27, as 
indicated by the observational results from SN la [1] and CMB anisotropics [2], and from 
studies of the large scale structure as well [3]. The origin of the cosmic dark energy is now a 
big challenge to astronomy and physics. There have been several types of models proposed to 
interpret the dark energy. The simplest one is the cosmological constant A. Since the matter 
density evolves as p m oc a(t)~ 3 and p\ ~ constant throughout the history of the Universe, the 
initial value of Pa/ p m has to be chosen with great precision to achieve p m ~ 0.37p A today. 
This is the coincidence problem [4]. Moreover, if A were interpreted to arise from the vacuum 
fluctuations of quantum fields [5], then there is a "fine-tuning" difficulty, i.e., why at present 
the vacuum energy of a scale 10~ 3 ev is so tiny compared to the typical scales in particle 
physics. Another model is the Steady State Universe [6], in which the accelerating expansion 
is driven by some C-field with negative energy, difficult to accept as a physical field. The 
effective gravity [7] is a model outside the scope of General Relativity, which has a task to to 
explain all observed features of gravity. The Born-lnfeld quantum condensate [8] is a model, 
which uses gauge fields as a candidate for the dark energy. One interesting speculations is 
that the dark energy may result from some dynamic field with a tracker behavior, i.e., the 
field is subdominant during early stages of expansion, later becomes dominant as the dark 
energy. Among this scalar kind of models are the quintessence [9], k-essence [10], phantom 
[11], quintom[12], etc. For a dynamic field model to solve the coincidence problem, it is 
required not to spoil the standard Big Bang cosmology, i.e., the Big Bang nucleosynthesis 
and the recombination must occur as usual, and the matter era must be long enough for 
structure formation. Moreover, the solution needs to be a stable tracker, insensitive to the 
initial conditions. But so far the scalar models have difficulty to fulfil these criteria [13]. 

In our previous work [14] [15] [16] a dynamic model is proposed, in which the effective 
YMC, up to 1-loop quantum corrections, serves as the dark energy. For quite generic initial 
conditions the model always has the desired tracking behavior, naturally solves the coincidence 
problem, in the sense that it satisfies those criteria mentioned above. When YMC is coupled 
with matter, EoS of YMC crosses -1 and takes w ~ —1.1, as indicated by the fittings of 
recent preliminary observational data [17] [18] [19] [20]. One may ask, does the model still 
work for more quantum corrections? As will be seen in this paper, the YMC dark energy model 
including the 2-loop quantum corrections has a tracking solution that solves the coincidence 
problem, like the 1-loop model. Through out the paper we use the units with c — h — 1. 



2. The Effective Lagrangian of 2-Loop YM Condensate 



The quantum effective YM condensate is described by the following Lagrangian up to 
2-loops [21, 22] 



L eff - -F 
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where, for the gauge group SU(N) without fermions, b = representing the 1-loop part, 
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and 7] = ~ 0.84 with &i = oV^C representing the 2-loop contribution, the parameter /t is 



the renormalization scale with dimension of squared mass, and F = — \F a ilv F aiiV = E 2 — B 2 
plays the role of the order parameter of the YM condensate. We consider only the case 
of "electric" condensate with F = E 2 . Here the dimensionless constant 5 is a parameter 
representing higher order corrections. Eq.(l) is equivalent to having the running coupling 
constant [23] 

with r = In |F/eK 2 | for L e ff = F/2g 2 (F). L e ff in Eq.(l) maintains the gauge invariance and 
the Lorentz invariance, and yields the correct trace anomaly [24] [21]. Eq.(l) is reminiscent 
of of the Coleman-Weinberg effective potential of scalar fields [25]. At strong fields, the 
corrections to the brackets in Eq.(l) are of order unity and, to the extent that they slowly 
vary, can be either absorbed into the scale mass k 1 / 2 , or represented by the parameter 5. As is 
explained in Ref.[22], renormalization group estimates are formally valid whenever the running 
coupling g 2 (F) is small in magnitude, (i.e., F/k 2 3> 1, giving g 2 (F) positive, and F/k 2 <^ 1, 
giving g 2 (F) negative). This suggests that there is a second weak-field asymptotically free 
regime, where Eqs.(l) and (2) give the leading behavior, and where g 2 {F) is negative, a feature 
of the effective YM condensate [14] [26]. By the way, the YMC introduced here is not the 
gluon fields in QCD, nor the gauge boson fields in the electro-weak unification. 
The energy density and the pressure, p y and p y , of YMC are given by 

Py = b -Fi [ T + 2 + r ] (\n\r + 5\ + -^)Y (3) 

^ = ^{r-2 + r ? (ln|r + 5|--A_)}. (4) 

So the form of the stress tensor T^ u of YM fields is consistent with homogeneity and isotropy 
of the Universe. The EoS for the YMC is given by 



When one sets rj = in the above expressions, the 1-loop model is recovered [15]. For 
comparison, we plot p y , p y , and w in Fig.l for both the 2-loop and the 1-loop order, respectively, 
where the variable is y = r + 1 = ln|F//t 2 |. We see that p y , p y , and w in the 2-loop order 
have similar shapes to, and higher magnitude than, those in the 1-loop one, respectively. At 
high energies y — > oo, p y and p y are positive, and the EoS of YMC approaches to that of a 
radiation, w — > 1/3, as is expected for an effective theory. At low energies y < 2.3, p y and 
w become negative, and at y < —0.8 the weak energy condition [14] [27] [28] is violated, 
Py + Py < 0, and w crosses —1. In particular, within the range of y ~ (—1.04, 54) attainable 
in our model of cosmic dynamic evolution, these physical quantities are all smooth functions. 
Although g 2 formally encounters a divergence at y ~ 0.4, nevertheless, the physical quantities, 
Py, p y and w, actually involved in our dynamical model are smooth, a similar situation as in 
the 1-loop model. The scale k can be fixed by requiring p y in Eq.(3) be equal to the dark 
energy density ~ 0.73p c , where p c is the critical density, yielding k 1 / 2 ~ 7.7 h^ 2 x 10~ 3 eV, 
where h is the Hubble parameter. At the moment we do not have an answer to the question 
why k is so small, so the "fine-tuning" problem is also present in our model. 

3. The Cosmic Evolution By the YM Condensate 

The spacetime of the universe is described by a spatially flat Robertson-Walker metric 



and is filled with the dark energy, the matter, including both baryons and dark matter, and 
the radiation. In our model, the dark energy component is represented by the YMC. To 
be concordant with the isotropy of the Universe, one may take an SU(2) YM field with a 
highly symmetric configuration Aq = and A\ = <j>(t)5 ai , where 4>(t) is a scalar function and 
a = 1,2,3 [29]. The overall cosmic expansion is determined by the Friedmann equations 



where p m and p r are the energy density of the matter and radiation components respectively. 
The dynamical evolutions of the three components are 



ds 2 = dt 2 — a 2 {t)8ijdx % dx^ 
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Py + 3-(Py+Py) = ~Yp, 
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Pr + 3-( Pr +p r ) =0, (10) 
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where T is the decay rate of the YMC into matter, a parameter of the model. Since T > 0, 
so the term Yp y in Eqs.(8) and (9) represents the rate of energy transfer from the YMC 
component to the matter one. Here for simplicity we do not consider the coupling between 
the YMC and the radiation. The sum of Eqs.(8), (9), and (10) guarantees that the total 
energy of the three components is still conserved. In terms of the variable N = In a(t), and of 
the re-scaled functions y, x = p m /\bK 2 , and r = p r /\bn 2 , the set of equations, (7) through 
(10), take the following form: 

N 2 = H 2 h 2 , (11) 
dy _ -% + 47y(ln \y-l + 5\ + ^) - j[y+l + yQn \y - 1 + 6\ + pf^)] 



dN y + 2 + r ] (\n\y-l + 5\ + ^_- w ^ w ) 



dx 7 
dN = ~ 3X+ h 



2 

y + l+r](\n\y-l + 5\ + 



y-1 + 6' 



(12) 
(13) 



^ = -4r, (14) 



where 7 = r/if, if = ^AnGbK 2 /3, h = ^x + r + ev\y + 1 + 77(111 |y - 1 + 5\ + pf+j)]. 
Looking at the forms, these equations are much more complicated than those in the 1-loop 
model with 77 = 0. Once the parameters 7 and S, as well as the initial conditions, are specified, 
the solution of this set of equations follows straightforwardly. 

As an example, we take the parameter 5 = 3 and the decay rate 7 = 0.5. The initial 
conditions for the set of Eqs.(ll) through (14) are chosen at a very high redshift Zi = 10 8 
during the radiation era. To ensure the equality of radiation-matter occurring at a redshift 
z = 3454 [2], the initial radiation and matter are taken as 

Xi = 3.68 x 10 22 , n = 7.25 x 10 26 . (15) 

To ensure the Big-Bang nucleosynthesis [30], the initial YMC fraction should be ~ 10% or 
less [16]. For concreteness we take it to be 

Vi < 54, i.e., ^ < 2.28 x 10~ 2 . (16) 

Pri 

The solutions for yi = 54 and for yi = 15 (i.e. — = 8.3 x 10~ 20 ) are explicitly shown in 
Figs. 2, 3 and 4. The dynamical evolution has a desired sequence of tracking solutions, i.e., 
during the radiation era the YMC follows the radiation as p y oc p r oc a(t) -4 , and then during 
the matter era it follows the matter approximately as p y oc p m oc a(t) -3 , and rather later 



around z ~ 0.6 it becomes dominant. Therefore, the matter era is long enough for formation 
of large scale structure. The fractional densities Q y = 0.73 and Vt m = 0.27 are achieved at 
z = (y = —0.97) as the current status of the Universe. Notice that this tracking behavior 
is always achieved for any initial p y i in Eq.(16) whose range stretches over ~ 25 order in 
magnitude. For a smaller initial p yi , say yi ~ 15, p y (t) only tracks p r (t) for a shorter period 
correspondingly, and then becomes the constant. Moreover, the dynamical equation for (y,x) 
at t — > oo has a fixed point (yf,Xf) = (—1.043, 0.114), which has been examined to be stable, 
i.e., the solution is an attractor. So, by the criteria mentioned in introduction, the coincidence 
problem is solved also in this 2-loop model. Furthermore, the EoS of YMC crosses w — — 1 
around z ~ 1.7 for the initial yi = 54, but for smaller yi the crossing occurs earlier. For any 
initial YMC, w ~ —1.14 at z = is obtained, which qualitatively agrees with the fittings of 
preliminary observational data [19] [20] [31] [32]. These behaviors are quite similar to those 
of the 1-loop model [15], even though the expressions involved in the 2-loop model are much 
lengthier. 

Parallel to the above of a constant 7, the desired evolution is also realized in an extended 
range of parameters, 7 ~ (0 ~ 0.6), and 5 ~ (3 ~ 7). By calculation, we find that larger 5 
and 7 yield a lower w at z — 0. For instances, at a fixed 7 = 0.5, 5 = 7 w ~ —1.18, more 
interestingly, at a fixed 5 = 3, 7 = 0.1 =>- w ~ —1.02, consistent with the A-model, indicated 
by WMAP [2], SNLS [17] and ESSENCE [18], and the crossing of -1 occurs rather late at 
z ~ 0.35, quite close to the fitting of Gold+HST sample [19]. The desired evolution can be 
also realized for 7 being a generic function of the YMC F. We have examined a variety of 
couplings, say 7 = 0.13e y//4 (2 + y)/(l + y) describing the YMC decaying into fermions and 
gauge bosons [33], and 7 = l.l%e v [15], etc.. In the simple case of non-coupling, 7 = 0, the 
dynamic evolution is similar to the coupling cases, but w does not cross, only approaches to 
— 1. Thus, if future observations do confirm EoS w < —1, then the coupling 7 > is needed 
in our model; otherwise, 7 = 0. 

4. Summary 

The dark energy model based on the 2-loop quantum effective YMC, under generic initial 
conditions and various forms of coupling, has a stable tracking solution of dynamics, thus 
naturally solves the coincidence problem. Moreover, with coupling to the matter, EoS of YMC 
crosses —1 and w ~ —1.1 at z = 0. In this model the matter era ends rather late at z ~ 0.6, 
so the matter era is long enough for the structure formation. Besides, since the YMC is 



subdominant during the early stages, so the nuclear synthesis and the recombination occur 
as in the standard Big Bang cosmology. These conclusions are the same as for the 1-loop 
model. With the Lagrangian being fixed by quantum corrections, the model depends on only 
two physical parameters: the coupling r ~ (0 ~ 0.6)H and the energy scale k 1 ^ 2 . As said 
earlier, since the "fine tuning" problem is not solved by the model, so the energy scale needs 
to be viewed as a new scale of physics. 
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Figure 1: p y , p y , and w as functions of the variable y = \n(F/n 2 ). The 2-loop model 
(solid) is compared with the 1-loop one (dotted) [15]. 



Figure 2: The evolution of energy densities in the 2-loop model with 7 = 0.5 and 5 = 3.0. 
Starting from Zi = 10 s the evolutions for two cases yi = 15 and 54 are given. For a 
whole wide range of initial YMC p y i = (10~ 27 ~ 10 _2 )p r j, there always exists a tracking 
solution, and p y {t) becomes dominant around z ~ 0.6. Due to the coupling, p m {t) will 
level off in future. 



Figure 3: The evolution of fractional energy densities in the same model as Fig. .2. At 
z ~ 3454 the equality of radiation-matter Q m = Q r occurred, and the values Q y = 0.73 
and Vt m = 0.27 are arrived at z — 0. 



Figure 4: The evolution of EoS in the same model as Fig. 2 and Fig. 3. In early stage at 
z — > 00, the YMC behaves like a radiation, w — > 1/3. Later on w reduces smoothly. Due 
to the decay of YMC, EoS crosses —1, and w ~ —1.14 at z — 0. 
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